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Abstract. Large differences between the properties of the known sample of cataclysmic variable 
stars (CVs) and the predictions of the theory of binary star evolution have long been recognised. 
However, because all existing CV samples suffer from strong selection effects, observational bias 
must be considered before it is possible to tell whether there is an inconsistency, which would imply 
a failure of the evolutionary model. We have modelled common selection effects and illustrate their 
^vq ' influence on observed CV samples. 
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INTRODUCTION 

O 

\Q \ Cataclysmic variables (CVs) are semi-detached binary stars with orbital periods {P or b) 

typically of the order of hours, consisting of a white dwarf primary accreting from a 
43 companion which is usually a late-type, approximately main-sequence star. See ll24ll for 

a comprehensive review of the subject. 

Binary star evolution is driven by loss of orbital angular momentum through gravita- 
tional quadrupole radiation, and, in at least some cases, the much more efficient mecha- 
nism of magnetic braking 0,111 III, EH. Observational estimates of the mass transfer rate 
(M) are on average much higher for long-period systems (P or b ^ 3 h) than for systems 
at P or b^ h J9[], suggesting that mass transfer is driven mainly by magnetic braking in 
long-period systems, while gravitational radiation dominates at P or j, < 3 h. 



The orbital period distribution of CVs is a useful indicator in the study of their 
evolution because loss of angular momentum, and the resulting mass transfer, changes 
P or b- This distribution is shown in the left hand panel of Fig. [T] Its most striking 
features are the period gap and period minimum: a pronounced drop in the number of 
systems at 2.2 h<P or b<2.& h and a sharp cut-off at about 80 min. T he p eriod gap 
is usually explained as resulting from disrupted magnetic braking fnl [Til The 
period minimum is a consequence of the response of a secondary of very small mass to 
continuing mass loss I12I1 . When the thermal time-scale of the secondary exceeds 
the mass-transfer time-scale, the star is not able to shrink fast enough (or in fact expands) 
in response to mass loss, so that the orbital evolution moves back through longer periods. 
CVs in this phase of evolution are referred to as 'period bouncers'. 

The rudiments of the theory of CV evolution have long been accepted, yet there 
are still significant discrepancies between the predictions of theory and the properties 
of the observed CV sample. Population synthesis studies predict that ~ 99% of the 
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FIGURE 1. The orbital period distribution of non-magnetic, hydrogen-rich CVs (left hand panel) | \ L 
Right hand panel: the same distribution divided into two magnitude bins, V > 16.5 (black) and V < 16.5 
(dark grey). The period gap is indicated by shading of the range 2.2h < P or j, < 2.8h in both panels. 



intrinsic galactic CV population should be short-period systems, and that about 70% 
should be period bouncers JEH. A glance at Fig. □ is enough to show that the first 
of these predictions is not reflected in the known sample of CVs. Furthermore, only 
a handful of known systems are likely period bouncers ificll . Clearly, if the standard 
theory of CV evolution is even approximately correct, the observed CV sample is not 
representative of the intrinsic population. However, since different types of CVs differ in 
intrinsic brightness and all surveys are flux limited, observed samples are not expected 
to reflect the intrinsic population. Apparent brightness and large amplitude variability 
are two factors that the discovery probability obviously depends on; since both intrinsic 
brightness and the frequency of large amplitude brightness variations decrease with M, 
short-period, low-M systems are necessarily under-represented in the known sample. 
The question is only whether these and other selection effects can account for the size of 
the discrepancy between the period distributions of the predicted and observed samples. 
The right-hand panel of Fig.[l]goes some way towards illustrating the importance of the 
apparent magnitudes of known CVs. 

We use a Monte Carlo approach to study the influence of simple but important 
selection effects on observational CV samples. 



THE COMPUTATIONAL METHOD 

A population of CVs drawn from the distribution in P or ^, (M), M\, and A/2 (where 
(M) is the long-term average mass transfer rate, and Mi and M2 are the primary and 
secondary mass) resulting from model pm5 of JfJ is distributed in a model galactic disc 
with exponential density profiles. The vertical scale heights are 120 pc for long- P n rh 
systems, 260 pc for normal short-P or ^ systems, and 450 pc for period bouncers ffllol. 
We take the galactic centre distance as 7620pc ||2j], and compute interstellar extinction 
by integrating the density of the interstellar medium [16] along each line of sight. 

In calculating the spectral energy distribution of each system, we account for the 
accretion disc, bright spot, white dwarf, and secondary star. The bright spot is modelled 



1 1 1 1 I 




log(P orb /h) 



FIGURE 2. Period distributions of magnitude limited samples. The magnitude limit is given in each 
panel, together with the fraction of the sample with periods above the period gap. Shading of the range 
2.63 h < P or i, < 3.35 h indicates the model period gap. The contribution of period bouncers is shaded in 
dark grey; the fine black histograms show DN that were simulated in outburst. 



as a blackbody. The flux contributions of the secondary and primary are found using the 
spectral type -P or b relation of [ 18] and the T e ff-(M) relation of lElll . We used the code 
of 1E2I1 to compute a 4D grid of accretion disc models over the variables (the rate 
at which mass flows through the disc), M\, the orbital inclination i, and the outer disc 
radius r c j. Most CVs have discs that are subject to a thermal instability; these systems are 
the dwarf novae (DN). The mass accretion rate through the disc in outburst (M^o) and 
quiescence (Mjq) are related by (M) = CM c /q + ( 1 — C)M^g, where C is the outburst duty 
cycle. Empirical estimates of the absolute V magnitude at maximum and a bolometric 
correction of — 1.8 yield the disc luminosity (L^); we set = 2R\Ld/GM\. ThenM^g 
follows from assuming C, and the probability of catching a DN in outburst in a single- 
epoch survey equals C. rj is set to 0.1 Ri for nova-like variables and DN in outburst, and 
0.5Rl for quiescent systems [3], where Ri is the Roche-lobe radius of the primary. 

SELECTION EFFECTS 

The single characteristic through which the most CVs have been discovered is large am- 
plitude variability (DN outbursts and nova eruptions). Most other CVs were found as 
objects with blue optical colours or X-ray emission, mainly in severely flux-limited sur- 
veys. A quantitative comparison between theory and observations requires a reasonably 
large observed sample with well defined selection criteria. We limit ourselves here to a 
mostly qualitative investigation of a few important selection effects. 



Optical flux limits 

Bias towards finding apparently bright objects is the most simple and important selection 
effect (see e.g. tl3J). Fig. El displays period histograms of samples with three different 
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FIGURE 3. Period histograms of systems with 
V < 20 (dotted), and systems with V < 20 and the 
U — B indicated in each panel (black). Dark grey 
histograms show the samples obtained when the 
same blue cuts are applied to the intrinsic colours. 
Almost all systems with V < 20 also have both 
(U - B) < -0.7 and U - B < -0.7, and most 
long-period CVs have (17 - B)o < - 1 . 1. 
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FIGURE 4. Period histograms of samples with 
the same U — B and (U — B)o cuts as in Fig. 3, but 
here applied to a volume limited sample. Dotted 
histograms show the complete sample; black are 
samples with observed colours satisfying the blue 
cut indicated in each panel, and grey shows sys- 
tems with intrinsic colours blue enough to satisfy 
the selection cuts. 



magnitude limits (V < 20, V < 16, and V < 14). In this figure (and all other period 
histograms we will show) the number of systems in each bin is scaled to reproduce 
a local space density of 5 x 10~ 5 pc~ 3 . The shape of the period distribution depends 
strongly on the magnitude limit. As is expected from the fact that long-period CVs are 
intrinsically brighter, the fraction of the total sample made up by long-period systems 
increases at brighter magnitude limits; it is 4.7%, 19%, and 33% for V < 20, 16, and 14 
respectively. 

Of course the known CV sample is not complete to 20th or even 16th magnitude. In 
fact, comparing the magnitude distribution of known CVs to a complete model sample 
indicates that the known sample is probably incomplete even at V = 13. Therefore the 
distributions in Fig. 1 should not be compared directly to the known sample. 



Blue optical colours 

Most known CVs have (U — B)q < —0.5 [1]. This is because high-M systems have high 
continuum colour temperatures, while low-M systems have the Balmer discontinuity in 
emission. However, in the systems with lowest M, the cool white dwarf is expected to 
dominate the optical flux. Blue selection cuts may therefore exclude the faintest CVs. 

Fig. 3 shows the effect of a magnitude limit together with different selection criteria 
in U — B. In every panel the dotted histogram shows all systems with V < 20, and black 
shows systems that also satisfy the blue cut. Only at U — B < —0.9 does the blue selection 
start removing a large fraction of systems from the magnitude limited sample. Fig. 4 
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FIGURE 5. The orbital period distributions of high galactic latitude (\b\ > 30°; black) and galactic 
plane (\b\ < 5°; dark grey) samples with magnitude limits V < 20 (left) and V < 16 (right). 



shows P or t, distributions resulting from blue cuts imposed on a volume limited sample 
(produced in such a way that all CVs have the same spacial distribution, so that this 
sample is representative of the intrinsic population). It is seen that even U — B < —0.5 
excludes systems near the period minimum; U — B < —0.7 introduces a severe bias 
against short-period systems. In the samples shown in Fig. 3, most of these systems were 
already excluded by the magnitude limit. A survey that selects objects for U — B < — 0.7 
is thus expected to be seriously biased against short-period, low-M CVs, but the colour 
cut introduces hardly any additional bias in a survey that is also severely flux limited, as 
is the case for existing t/V-excess surveys. 

The grey histograms in Fig. 3 and 4 display samples with the same blue cuts as 
shown by black histograms, but here the cuts are applied to the intrinsic colours to 
illustrate the effect of reddening. Notice that practically all long-period systems have 
(U — B)q < —0.9, but some get excluded from observed samples because they are on 
average more distant than short-period systems, and thus reddened. This is amplified by 
the fact that they are also more concentrated towards the galactic plane. 



Restrictions in galactic latitude 

Surveys at high galactic latitude hold some promise of uncovering samples that are 
volume limited. In a galactic plane survey, on the other hand, intrinsically bright objects 
at large distances will be included in disproportion to their space density. The direction 
of this effect is obvious; we simply show its magnitude in Fig.|5l with period histograms 
of high (black) and low (dark grey) galactic latitude samples with two magnitude limits. 
As expected, the ratio of short- to long-period systems is smaller in the galactic plane 
samples, and this effect is larger for the deeper samples. 

In a survey covering the area \b\ > 20°, less than 10.0% of CVs in the survey volume 
are detected above a magnitude limit of B < 20. This implies that even the CV sample 
produced by a deep, high-Z? survey such as the Sloan Digital Sky Survey Il20ll will by no 
means be volume limited. Note also that, since the vertical galactic distribution of CVs 
differ for sub-populations with different typical ages, a volume limited, high-Z? sample 
is in any case not the same as the intrinsic CV population. 



CONCLUSIONS 



We have presented a Monte-Carlo technique to model selection effects in observed sam- 
ples of CVs. Given a theoretical intrinsic population and any survey with well defined 
selection criteria, this technique can be applied to predict the observed population. 

For practically achievable flux limits, no flux limited sample is representative of 
the intrinsic CV population; it is therefore quite inappropriate to compare an observed 
sample with the predicted intrinsic population. The properties of the period distribution 
of a magnitude limited sample depend strongly on the magnitude limit. A selection cut of 
U — B < — 0.7 introduces a serious bias against detecting short-period systems, but this is 
unimportant in the presence of a bright flux limit. Finally, magnitude limited surveys at 
high galactic latitude are expected to produce samples very different from galactic plane 
surveys, because they detect a larger fraction of all systems inside the volume defined 
by the galactic latitude range. 
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